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Abstract 

The experimental demonstration that neutrons can reside in gravitational 
qimntum stationary states formed in the gravitational field of the Earth 
indicates a need to examine in more detail the general theoretical 
properties of gravitational eigenstates. Despite the almost universal study 
of quantum theory applied to atomic and molecular states very little work 
has been done to investigate the properties of the hypothetical stationary 
states that should exist in similar tj^es of gravitational central potential 
wells, particularly those with large quantum numbers. In this first of a 
series of papers, we attempt to address this shortfall by developing 
analytic, non-integral expressions for the electromagnetic dipole state -to- 
state transition rates of charged particles for any given initial and final 
gravitational quantum states. The expressions are non-relativistic and 
hence valid provided the eigenstate wavefunctions do not extend 
significantly into regions of strong gravity. The formulae may be used to 
obtain tractable approximations to the transition rates that can be used to 
give general trends associated with certain types of transitions. 
Surprisingly, we find that some of the high angular momentum eigenstates 
have extremely long lifetimes and a resulting stability that belies the 
multitude of channels available for state decay. 

PACS numbers: 03.65.Ge, 03.67.Lx, 03.65.Db, 04.60.-m, 95.35.+d, 
04.90.+e 



1. Introduction 

The bound stationary states of electrons in atoms have been studied extensively 
for many years but little has been done to investigate the properties of neutral or 
charged particles occupying the hypothetical stationary states in gravitational 
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potential wells. This is important because there are no strong scientific grounds 
to deny the existence of gravitationally bound quantum eigenstates (including 
macroscopic ones) and also because recent experimental work by Nesvizhevsky 
et al [1,2] has demonstrated the physical reality of such gravitational quantum 
stationary states in the Earth's gravitational field. It is interesting therefore to 
speculate on whether it might be theoretically possible for relatively pure 
stationary gravitational eigenstates to exist naturally elsewhere in the universe. If 
such states were to exist, it would clearly be important to have theoretical 
information on their expected properties. 

This paper aims to set up a simple quantum model using a central gravitational 

point potential, derive its set of hypothetical gravitational eigenstates and begin 
an analysis of the properties and interactions of some of these states. Specifically 
we develop formulae for the form of the dipole-allowed state-to-state transition 
rates in the general case of charged particles occupying stationary states with 
arbitrarily large quantum numbers. Although mathematically the eigenstates of 
such a system are theoretically analogous to those of an equivalently bound 
electrical system such as the hydrogen atom, differences in the sizes of the 
quantum numbers and scales involved, combined with the unusual properties of 
the eigenfunctions themselves, result in some of the gravitational eigenstates 
having unusual and unexpected properties. Deriving these properties and 
understanding reasons for the differences between the atomic and gravitational 
cases is often not trivial. The long term aim is to determine under what 
conditions long-lived, relatively stable, gravitational eigenstates can exist, and 
further, to make predictions about the properties of possible (necessarily large 
scale) hypothetical structures composed of an ensemble of such states. If such 
relatively pure gravitational eigenstates do exist in the universe then it is 
possible that they could have significant consequences for Astronomy [3, 4]. 

The simplest realistic gravitational system consists of a symmetric potential 
produced by a sufficiently large mass M so as to have well bound eigenstates 

occupied by particles of atomic mass scales. We therefore consider the 
gravitational quantization of a small particle mass (« Af) in the point 

potential field provided by the large, electrically neutral, mass M. To simplify 
matters fiirther we consider the states whose radial extents are limited and have 
appreciable amplitude only over regions which are sufficiently distant fi^om the 
centre of M so that the gravitational field is weak enough to ignore relativistic 
effects. Under these conditions the Schrodinger equation is trivially analogous to 
that for the hydrogen atom and may be written as 

Gm M Qy, 

-—V 11/ i^ii/ = iTl— (1) 

2fj. r ot 

where /j. is the reduced mass and the other symbols have their normal meanings. 
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The solutions are likewise analogous to those of the hydrogen atom and may be 
immediately written down, the eigenvalues being 



En=- 



(2) 



With the introduction of the parameter - 
eigenfunctions u„(r,t) are 



GjuntpM 



the corresponding 



Wn,/,m(r,0 = ^„,/('')j/,m(^>^) 



(3) 



where 7;^(^,^) are the normalized spherical harmonics and 
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In equation (4) N„i 
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In (« + /)! 



is a normalizing constant and 
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(5) 



are the generalized Laguerre polynomials in their standard form. 

In most traditional situations where a quantum approach would be anticipated, 
the resulting eigenstate binding energies are unrealistically small because the 
gravitational force is relatively very weak. For example the binding energies E„ 
of two neutrons (ignoring spin) is much smaller (e.g. Ei ~ 10"*' eV) than the 
magnitude of typical random background field fluctuations (e.g. E ~ 10''^ eV for 
cosmic microwave backgrovmd radiation). The binding energy increases with 
larger masses, but for small quantum numbers the physical extent of the 
quantum probability density distribution can become much smaller than the 
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corresponding quantum probability distributions of the individual isolated 
masses themselves making up the structure (M~ m ~ 10"'"' kg implies an « = 1 
eigenstate 'size' of ~ 10"'^ m compared to a typical 'size' for the individual 10"'^ 
kg masses of -10"^ m for a mass density ~1 kg m'^). The two particle gravity- 
only Schrodinger equation is inapplicable in this situation since the interactive 
effects of the electron clouds in the two masses would then dominate. We 
therefore consider the theoretical description and properties of very large, 
possibly macroscopic, stationary, well bound states that correspond to large 
quantum numbers in the (relatively) weak regions of deep gravitational wells. 

The notion of large quantum states and macroscopic quantum phenomena is not 

new. Bose-Einstein condensates. Superconducting Quantum Interference 
Devices, and the many experiments involving quantum connectedness all 
demonstrate existence of macroscopic quantum phenomena [5-8]. Of course 
there are understandable reasons for certain quantum phenomena, such as 
interference effects, not occurring on large scales. In the case of quantum 
interference, for example, background electromagnetic or gravitational field 
fluctuations WiW generally produce interactions that result in decoherence on 
macroscopic scales. However these only introduce randomly phase shifts that 
occur independently to the individual members of a quantum connected pair of 
virtual states [9] thereby destroying the quantum coherence and the subsequent 
interference. There is no experimental evidence that this type of macroscopic 
decoherence should result in the indiscriminate rejection of all macroscopic 
quantum phenomenon. 

Not only will we be concerned with high-M and high-« values, but also the 
relatively high angular momentum states (high-/ relative to n) because for these 
states the radial extent of the probability density 'shell' can be relatively limited 
(and conversely, low-/ relative to n implies large radial eigenfunction spread). 
When the central mass is small, then, as the n value increases, the binding 
energy becomes too weak for effective binding. This occurs long before the 
effective position of a high-/ eigenstate is able to outstrip the physical size of the 
central mass M in question, but can be easily achieved if the central mass is large 
enough. The lowest central mass limit at which this phenomenon becomes 
feasible depends on the assumed density of the central mass and on the critical 
value that one adopts as an 'effective' binding energy. However, for realistic 
mass densities of 10° - 10^ kg m"', and a lower limit binding energy of around 1 
eV (i.e. E„ » average background CMB photon energy) the required central 
mass M is of the order of 1 0^" kg. This explains the difficulty in producing or 
observing gravitational eigenstates on laboratory or microscopic scales. It is of 
interest to note that the binding energies of the neutrons observed in the work of 
Nesvizhevsky et al (quantized in an asymmetric wedge potential of a 'hard' 

25 

mirror and the Earth's ( M < 1 kg) gravitational field rather than a central 
spherically symmetric one) had binding energies of only 1 .4 to 4. 1 peV for « = 1 
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to 4 yet still demonstrated relatively large 'sizes' by traditional quantum 
standards (-0.05 mm for « = 4) [2]. 

The solutions to (1) are not trivially dealt with by analogy with those of the 

hydrogenic wavefunctions because of the necessarily large quantum numbers 
involved in the formation of energetically realistic gravitational eigenstates. Any 
interaction properties that depend on the relative energy level spacing and the 
overlap integrals of states involved can be drastically different from those 
properties resulting from low-« quantum state overlap integrals. Large quantum 
number wavefunctions of electrically bound electrons in hydrogen need to be 
rarely, if ever, considered in the atomic case because their binding energies are 
too small to remain stable and bound for any reasonable length of time. 
Describing and dealing with states having high-«, /, m values is therefore a 
primary goal of the present work. When n is large the complexity of the 
individual states and the number of available eigenstate transitions makes the 
general study of the intrinsic properties and interaction rates of any large array of 
gravitational eigenstates (referred to from now on as an 'eigenstructure') 
exfremely difficult. 

We inttoduce a schematic (figure 1) showing the fraditional quantum parameters 
n, I and m. As with the atomic case, the total angular momentum parameter / and 
its z-projection sub-levels m, range from to n — \ and -/ to / respectively. 
Vertical lines represent states of constant / and horizontal lines states of 
constant n . Each soUd circle in figure 1 represents the 2/ + 1 gravitational z- 
projection sublevels. It is convenient to infroduce a quantity p(=n-l) that 
then ranges from 1 to n. The reason for this is that it will turn out that the states 
of interest here will have p values that are small relative to n and formulae for 
these eigenfunctions are more conveniently written in terms of p. The diagonal 
dotted lines of figure 1 represent lines of constant p and all states lying on one of 
these diagonal lines all share a common p value, beginning with p = l on the 
leftmost diagonal. The parameter p reflects the general spatial profile of the 
radial eigenfimction component: p = l corresponds to a single-peaked, radial 
eigenfimction component with no zeros, p = 2 a two-peaked, one-zero ftmction, 
p = 3 a three-peaked, two-zero function and so on. 
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l=n-2 l=n-4 



increasing / value 



/ = 




increasing 
n value 



Figure 1. Schematic representing the high n, 1 ,m -valued stationary 
states, drawn to emphasize the parameter p = n — l . Each solid circle 
on the diagram represents 2/ + 1 z -projection substates. 

To aid the physical understanding of the behaviour of large quantum-valued 
eigenstates we end this introduction by giving examples of the eigenfunction 

probability functions w„ (r,t)u^{r,t) of typical gravitational systems for 

instructive n, I, m values (essentially a revision of traditional hydrogenic 
systems) and note some relevant features. The spherical harmonic components 
7/ „, {d,<^) of the eigenfunctions are explicitly given by 
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(-1)'" \ 2l + \ (I -my. 
I'n V 4;r {l + m)\ 



(2^)!(cos(^)) 
k\{l-k)\ {2k-l-m)\ 



2k— I— m 



\m{0,<l>) 



(sin 6')' 



m 




(6) 



The azimuthal dependence on ^ involves only functions of the form e"*"""^ and 
as with the atomic case, gravitational eigenfunctions result in probability density 
functions that are completely delocalized in the (j) direction. Hence, one simple 
representation of the density functions of the stationary states presents them as 
three dimensional revolutions around the ^ = axis of two dimensional r, 
polar coordinate plots shown here by defining a radial-polar probability surface 
density associated with the probability of finding a particle between (r, rO) and 



element). Examples of the radial-polar probability density dependence on n, I, m 
may then be drawn (taking an arbitrary scale parameter Z?o^l) as in figure 2. 
They are shown here to emphasize that (a) the number of peaks in the density 
function in the radial and polar directions is given by p and (/ - w + 1) 
respectively and (b) the radial 'thickness' or extent of the state probability density 
reduces with decreasing p. We also note the other well known hydrogenic 
eigenfunction properties that the average radial position of the density function 
increases as n increases and also that, for a given normalized set of states with 
the same n, I values, the summed probability density over all possible m values 
gives a symmetric distribution whose integral over all space is 2/ + 1 . For arrays 
covering sufficiently large n and / values, the sum over a sparsely populated set 
of m states will also show spherical symmetry provided that the occupancy is 
randomly distributed. The degree of localization in the r and 6 directions varies 
with values of m and / but, predictably, none show the degree of locaUzation 
exhibited by traditional orbiting particles. It is clear that the quantum 
representation of any spatially-localized orbiting particle in terms of the 
eigenstate array will necessarily involve a summation over many eigenstates and 
that the tighter the localization constraint the more expansive is the required set 
of eigenfunctions. 



{r + dr,r{6+d6)) as 



d^ (where c/F = infinitesimal volume 



drdO dV 
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Figure 2, 

(a) Peaks in the radial-polar density profile reflect the values of p 
and l-m+\. bQ= \, n = A5,{jL) p = 15, /-m+1 = 6; (ii) p = 15, /-m+1 
= 21; (iii)j9 = 5, l-m+l =30; 

(b) Radial 'thickness' of density function decreases with 
decreasing p. 

Z?o = 1, « = 50, m = 3, (i) j3 = 10; (ii) p = 3; (iii) p= I; 

For large quantum numbers calculation of the various overlap integrals for 
spontaneous dipole decay is particularly difficult and it is useful to derive non- 
integral formulae for these which in turn lead to more manageable formulae for 
calculating state-to-state transition rates and ultimately state lifetimes. This is 
done in section 2. Despite the more manageable forms of the non-integral 
formulae for dipole decay, it becomes necessary to introduce approximations for 
specific cases. One case of relevance to the present work is that where one or 
more of the states has a p quantum value of 1 (see figure 1). Section 3 deals with 
the development of two useful approximate formulae which enable calculations 
for these types of transitions. Lastly we conclude with a discussion of the 
findings of this paper, their relevance to speculations on the existence of 
naturally occurring gravitational eigenstates in the universe, and the need for the 
development of further approximations that will be the subject of another 
companion paper. 



2. Exact non-integral state-to-state transition rate formulae 
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The interactive properties of particles in gravitational eigenstates will be 
determined by an overlap integral involving a relevant interaction potential. It is 
not possible in the present treatment to cover all possible types of interaction 
potentials. In the atomic case the most basic questions concern the stability and 
decay of the electronic states. We therefore restrict attention initially in this 
paper to the calculation of decay rates and the intrinsic stability of the states. 
Subsequently it will be shown in later papers that the mechanisms determining 
some of the important properties of certain of the states discussed here are also 
applicable to more general interactions because of intrinsic characteristics of the 
eigenstates themselves, and are not greatly dependent on the form of the 
interaction potential that is part of the overlap integral. State-to-state transition 
rates and state decay rates may therefore enable the estimation of rates for other 
types of interactions as well. 

The central mass is assumed neutral so that eigenstate radiative decay takes 
place via emission of either gravitational radiation or, if the eigenstate particle is 
charged, a combination of gravitational and electromagnetic radiation. It has 
been shown [4] that for the typical conditions considered here radiative decay 
through gravitational radiation is insignificant compared to that produced 
through electromagnetic dipole decay. We therefore derive formulae for the 
state-to-state electromagnetic dipole decay rates that would be applicable, for 
example, to the two most stable particles, the electron and proton in various 
eigenstate configurations. 

If the eigenstructure is sufficiently densely occupied, then fields of the eigenstate 
particles themselves will also contribute significantly to the potential as seen by 
any individual eigenstate particle. However in a regular eigenstructure array with 
equal numbers of alternating positive and negative charges, it may be shown 
that, from a classical perspective [4], the electromagnetic component of the 
potential felt by an individual charge will, for a large central mass and low 
particle densities like those encountered in some interesting astrophysical 
situations be negligible compared to the gravitational potential. For example, 
consider the charged eigenstate particle in (1) to have charge q and be 

embedded in a uniform array of other positive and negative charges qj (= ±q) , 

then the electrical contribution to the potential energy term from these other 
charges maybe estimated from the analogous solid-state formula 



other alternating charges, a is the Madelung constant and s is an average 
separation distance (analogous to a lattice constant). For typical present-day 
galactic halo conditions where the average occupied proton/electron particle 

eigenstate density is ~10^m"^ and the central mass ~10''^ kg, the above 




where r, is the distance to each of the 
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formula gives the electrostatic potential energy of a typical proton as ~ 10 
times smaller than its gravitational potential energy. Limits on the assumption of 
a neutral central mass M may be obtained by comparing the gravitational 
potential energy term in (1) with an electrical potential energy term when the 

central/enclosed mass M has a net charge Q , for a proton in a typical 

eigenstate. In this case these terms are GMm I r and Qq I Ans^r respectively. 

Substituting typical halo values of Af ~ lO'^^kg and r ~ 10^' m, shows that the 

central mass M would need to have a net charge g ~ 10^"* C for the electrical 

potential energy term to be of comparable magnitude to the gravitational energy 
term. The existence of such large charge displacements within the galaxy would 
seem unlikely, but this comparison provides a figure for determining a limit on 
the validity of the neutrality assumption in terms of the overall galactic charge 
displacements in the present treatment. 

In the present work therefore we ignore any electrical contribution from the 
charged eigenstate particles themselves to the overall potential as seen by any 
one particle and also take the central mass M as neutral, subject to the 
limitations on this as expressed above. Furthermore as regards to eigenstate 
particles' global contribution to the gravitational potential, since we will be 
looking predominantly at high n , / states where the thickness of their shell-like 
radial probability distributions is relatively small, the total effect on the potential 
of the other eigenstate particles may, in this case, be incorporated into an 
effective net central gravitational potential, analogous to that used in the 
Hartree-Fock approximation in atomic physics. The condition that this net 
central potential is weak over the physical extent of the eigenstate is given by the 

general condition for weak gravity s = GM I rc^ « 1 . For particles in the 

region of a typical galactic halo for example, M~10''^kg and r~10^'m, 

giving £ ~ 10"* , so that this condition is well-satisfied. Since halo density scales 

as 1 / , M scales directly with r and this condition is satisfied at lower r as 

well, until for r ~ lO^'^m, and typical central galactic black hole mass of 10^^ 

kg, £ — 10"^. Thus for the low-/? (<10'°), high-« (>10^') eigenstates 

considered here, whose radial extents are <10'''m, and whose radial positions 

begin at an r at least as large as r ~ 10"'(» lO''*) , the condition of weak 

gravity, the 1 / r approximation to the variation of potential, and the neutrality 
condition discussed above, are simultaneously satisfied. 



10 



Gravitational eigenstates in weak gravity I: . 



The state-to-state transition probability for radiative dipole decay A^J { = 

Einstein A coefficient)^ is analogous to the atomic case. For a transition «, Uf 
from initial state |/) to final state \f) this is [10] 



"^j ~ ~~ rr ^'^ 

where e is the electronic charge, Sq the electrical permittivity of free space, 
|^/|er|/^| =n,j the absolute value of the dipole matrix element for spontaneous 

decay for the transition / to f , C0ff = 

angular frequency corresponding to the transition Mo f , ju the reduced mass, 
and the other symbols have their normal meanings. For a state / with multiple 
decay channels, the reciprocal state hfetime 1 / is then the decay rate sum over 

all possible decay channels k: 1 / r,- = ^ = . 

k 

In calculating 1 / there are two problems: the number of transitions k that must 
be summed over may be extensive for some states, and it may be very difficult 
to calculate n,y(=n-^) for the particular |/) and | A:) involved in each of the 

individual (co^i^ may be calculated fi^om the eigenvalues in a straight 

forward manner however.) Explicitly Yl(f is written as: 



4 



(8) 



where FI;-^, Ilyy and n,-^ are the x , y and z components of the vector inside 

the modulus in (8). Dipole radiative decay occurs via transitions involving 
Aot = (implying ll^y^ = Il^yj, = 0) or tsm-±\ (implying n,y^=0) and 

A/ = +l (implying transitions must take place between adjacent /columns in 
Fig. 1. 



' Higher order decays show similar trends to that of dipole decay but rates are several orders of 
magnitude smaller for each corresponding unit increase in multi-pole order. 
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n,-^, n,-^ and n^-^ may be explicitly written as 



"rx = ^nf,if (^rRnUi Yii^Mr^ sin(6')cos(^)sin(6') J^6'(/r 

(9) 

%=Jo Jo Jo sin(6')sin(<z))sin(6')J<zic/6'afr 
= ^lo^hf ^nijAj^C ^V^rnf \mi sin(^) sin((^) sin(^)t/^^ 



= ^^R ^e^y 



= ^lo^nfjf ^nijAl^C ^IfM \rni COS(0)sm(0)d^0 
= ^lRhd,z 



(10) 



(11) 



where we have further split the integrals into their radial ( 
^«=Jo Kf,if('') Kiji(r)r dr) and angular ( 

h^x =10^^^'}'"^ \mi^^(^)^os((^>m(0)d0d<f> , etc.) components. 

Explicit formulae for the radial and angular integrals in (9), (10) and (1 1) may be 
then obtained using (4) and (6) respectively. 

The explicit forms of the angular integrals Ig^^ , Ig^y and Ig^^ vary depending 

on whether Am is +1, -1 or and whether A/ is +1 or -1. We briefly show here 
the technique for derivation of the case where the initial state is J/,^, = J;^ and 

final state isYij-^^j- =Yi_i„j, that is the combination Am = and A/ = -l, 

derivation of the other combinations following along similar lines. In this case 
7(51^^ and I^^y are zero and we need only calculate Ig^^ , 

C TU C '2.7Z -if, 

Jo Jo Assuming for the moment that / and m 

are both even, using Rodrigues' formula and the binomial theorem, and 



12 



Gravitational eigenstates in weak gravity I: . 



adjusting the summation limits appropriately, gives the initial state J^,^, = 1/ „ 
as 



4;r {l + m)\ 



(sin^re""<*x 



\{l-m)l2 
fc=0 



(2/-2^,)!(cosg/^-^^'-"') ' 
2'(/-2^,-m)!(/-^,)!^/! 



(12) 



and the final state i//^^/ = i/-i,m 



(m) 



(2/-l)(/-w-l)! . 



4;r (/ + m-l)! 



(sin 61)™ e' 



{l-m-\)ll 



. (2/-2A:^-2)!(cos^)^ 



(/-2fc^-m-l) 



(13) 



(l-2kf -m -1)1(1 -kf -1)1 kj ! 

where it is understood that the summation variables extend only to floor 
integer values for limits which become half integer. 



then becomes 

f 



Jo Jo 



An (/ + /77)! 



(2/-l)(/-»j-l)! 



\k (/ + W-1)! 
( 



(sin^)"'e- 



\{l-m)ll 

s 

/t,=0 



\(/-2fcj-m) 



_ (2/-2^,)!(cos^)^ 

l\l-lk^-m)\{l-k^)\k^\ 



{l-m-\)l2 
kj=0 



A,_(2/-2^^-2)!(cos^?) 



(/-2fey-W-l) 



T'\l -Ikf-m- !)!(/ - A:^ - 1)!A:^ ! 



cos ^ sin ^ 



V 



d0d<p (14) 



The product of the summations may be then combined into a double summation, 
the powers of sin^ and cos^ collected, the integrals brought inside the 
summation and the result simphfied to give 
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1 1(2/ + l)(2/-l)(/ -/;;)!(/ -/»-!)! 
(/ + w)!(/ + w-l)! 



(_l/-+^/(2/-2^,)! , 
2^^-'k,\kf\{l-k,)\' 







(l-m-\)l2 


(l-m)/2 

z 




V 



(2/-2A:^-2)! 



il-2k.-m-\y.il-2ki-my.{l-k.-iy. 



f 

{2l-2ki-2kj-2m) 



dO dcj) 



(15) 



The integral over ^ is just 2k and the integrals over 9 are of the form 
['^ {?,m9)^^^^ (co^O)^^ dO v^fith r and 5 are arbitrary positive integers. These may 

JO 

be evaluated using the recurrence relation 

(sin ef'^^ (cos ef' dO = (sin Of'^^ (cos Of''^ dO 

■"0 2r + 2^ + l-'0 

and the fact that 



f^/ • m2r+i ^ csc(;rr) ^/^^(l + r) v • * 

(sm ^) t/^ = - = -— ^ (for positive integer r 

JO r(-r)r(3< + r) r()/ + r) 

to give 



•) 







[ ' 1 


U'' + 25 + l j 


lv2r + 2^-l ]"' 





i-lr(sin^)-V^ 



2^^(25)!r!(2r + l)!(r + 5)! 
'()^ + r)(2r)!s!(2r + 2s + l)! 



(16) 



We set r = m and s = l-k^-kj■-m in (16), substitute these into each of the 



^2l-2k-2kf-2m 



expanded and collected product sum terms (sin ^)^'""^^ (cos ^)^' ^'^^ 
(15). 



of 



This yields a final expression for Iq^^ , so that the angular integrals for this case 
may be finally reduced to 
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h(l>z = 

f 



^im-iM^, (2/ + l)(2/-l)(/-m)!(/-m-l)! 







(/-OT-l)/2 

■ z 


(/-m)/2 

z 







(-1) ' ^ (2/ - 2^,.)! (2/ - 2^. - 2)!(/ - - A:.)! 
yt,- lyt^ !(/ - kj - !)!(/ - A:.)!(2/ - 2A:,. - 2A:^ + 1)! ' 

{2l-2m-2ki-2kf)\ 
{l-2kf-m-\)\il-2ki-m)\{l-ki-kf-m)\ 



(17) 



Similar alternative expressions can be obtained when the quantum numbers / 
and m are odd or combinations of odd and even, or when the final state is / + 1 
rather than Likewise expressions for /^^^ and Ig^y can be obtained for the 

cases where Aw = +1 and in this way explicit expressions for all dipole 
transition angular integral components may be obtained. 

The correctness and exactness of (17) can be verified by comparing cases where 

the values of / and m are small enough to also allow alternative calculation by 
direct integration. For example for / = 30 w = 6 , both integration and (17) yield 
the exact result of 12.^6/3599 . 



* 3 

The derivation of the radial integral, ^«-Jq ^nf,ifi'')^ni,iA^)^ is 

somewhat more involved but, unlike the angular integrals, only takes one of two 
possible forms depending on the angular momentum / + 1 or /-I of the final 
state. We present here the latter case taking the initial state as 
{n^, lj) = (nj, 1^ -n^ - p) where p -n^ as defined earlier, and final state as 
(nf,lf) = (nf,lr=li-l = ni-p-Y). The initial and final states written in 



terms of p then respectively become, using (4) and (5) 
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R 



n.J 



( 1 V 



{p-\)\{2ni-p)\ 



2«,- 



exp 



(-1)' 



f^,{p-k,-l)\{2n,-2p + k, + l)lk,l 



and 



R 



n y ,1-1 



z 



' 2 ^ ' 



2«/ 



fc^=o (nf-ni + p-kfy.i^2ni-2p + kf -l)\kf\ 



exp 



2r 



(18) 



I]f can be then written as 



Cx X 



p-1 

z 

fc=0 



V V 



(■00 



nj+rif 
"in/bo 



2ni-2p+k,+kf+2 



dr 



nS/'k,\kf\{p-k,-\)\{nf+p-ni-kf)\ 

1 

{2ni-2p + kf-\)\{2ni-2p + ki+\)\ 



J) 
(19) 



where 
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C- 



In,- 



Since 



poo 



dr 



^[2{ni-p) + k,+kf+2)\ 



both"/ 



/ poo ^ 3 I 

the expression for I = ^ (r) R^- /■ (r) r dr\ can therefore finally be 

written as 



Ir = 



((p-l)!(2«, + -p-\y)^ X 



(20) 



kf=0 



p-1 

z 

£=0 



(2«,.-2;? + ^,.+A:^ + 2)! 
(2«,-2;, + A:^-l)!(2«,-2;, + ^,+l)! 



V V 



Again the exactness of (20) may be checked by substituting low values of , 
nj- and p directly and comparing these with results obtained by direct 

integration of ^nfjf ('') ^niji ^ . For example, if n, = 6 , Uj- = 3 , p-4 
(so that = 2 and // = 1 ) both expressions give a value of 4096^/70 Z>q / 19683 . 



Combining (7), (17) and (20) then enables any state-to-state transition rate to be 
calculated for A/ = — 1 and Am = . As a check on the validity of these 



17 



Gravitational eigenstates in weak gravity I: 



expressions, they may be easily converted to their equivalent atomic 
counterparts. For example the decay rate for the =2,1. = \^n^ =i,lf =0 

transition in hydrogen may be calculated using the present treatment by 
substituting parameter values relevant to this atomic transition yielding 
6.2x10^ s~' compared to literature values [10, 11] of 6.3x10^ Decay of 
the "Rydberg" atomic hydrogen «, =10, /, =0 and =10, /, =1 states are 

more complicated because they involve addition of rates over a greater number 
of decay channels and also transitions involving other A/ and Am values ( 
A/ = ±l and Am = 0,±l). Using equations (7), (17) and (20) and other 

equivalent expressions the present treatment yields lifetimes for the hydrogen 
Rydberg states «, =10, /, =0 and =10, =lof 2.0 and 0.19 //s compared 

with the values obtained from [12] of 2.1 and 0.21 jus respectively. 

3. Dipole matrix element and transition rate approximations for transitions 
where at least one state has p = l 

The interesting states are those that have very low decay rates and that are 
weakly interacting as these would be the ones most likely to exist naturally. We 
are concerned therefore to know upper limits on the decay and interaction rates 
of particular types of transitions. In the case of dipole decay these depend on the 
values of H^y and ft>,y in (7). The value of is critical in determining the 

decay rate because it depends on r that can become very large. The smallness of 
angular components of Tl^j therefore becomes an important consideration in 

limiting the transition rate only if Ig is not small. Consequently it is the value of 
Ijf that we wish to primarily investigate here, because of its dramatic effect on 
Hff and on the resulting transition rates. 

No approximations have been made in the derivation of /jj in (16) and (20) so 
that these equations enable the possibility of obtaining tractable expressions for 
the radial component of the dipole matrix elements corresponding to any dipole 
transitions within the conditions set out earUer in this paper. For example 
substituting «y =n, -1 and p = \ into (20) gives an explicit expression for the 

radial component of transitions of the type B to A ( 
(«;• = 1) ^ («y = nj - 1,1 = Til - 2) ) shown in figure 1 as: 
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f 



|;i?;,^(r)r^i?„,,(r)^r=22"-^jf^^ ^(2«, -if (2«, -2) (21) 

which, for large n^ becomes ~ ft^n. . Again the agreement of the formulae for 
these low-/?-to-low-/? state transitions may be verified when the values of n, I and 
p are sufficiently low for to be manageable by direct integration of the 
polynomial representations. 

A similar dependence on Z>o«, exists for adjacent state to state transitions on the 
p — 2 diagonal and hence yields radiative decay rates essentially the same as 
that for transitions of the type B to A. It can be shown using (20) that whenever 
transitions take place between two '/j-tuming point' Laguerre eigenfunctions, 
that is transitions that take place along the same /j-diagonal so that «, -^n^-\ 

and — > - 1 , then the decay rate is h^n^ provided » p . 

Transitions like C to A of figure 1 involve overlap integrals where the radial 
component wave functions have very different shapes (A is a 1 -turning point 
function while C is a 2-tuming point function). As a result, it would be expected 
in this case that the radial part of the overlap integral would be much smaller 
than the B to A type transitions. This is indeed the case and (20) gives 



\lKf,lfiry''R^^^l^{r)dr= 



>2«,-2, 



(2«,-2f 



V(2«,-2)(2«,-3)(2«,-4) 



(22) 



which reduces to the approximate value of h^n^^ / %/2 again provided » p . 

A general formula for the value of 1^ may be obtained for a more general 
transition such as D to A of figure 1, which originates from an arbitrary p- 
tuming point radial Laguerre polynomial state D («,,/, =n--p^ and ends on a 
1 -turning point state like A («y =ni- p,lf =nj- - I , lying on the left diagonal). 

In general for such transitions it is possible for p to be much larger than 1 (but 
still much less than «j ). It may be shovm that in this case (20) reduces to a single 
summation 
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Kj-p 



({p-l)\{2ni-p)\(nf+p-ni)\(ni+nf-p-l)^Yx 

{2ni-2p + ki + 2)\[-2nff 



L=0 



(n,+nf)-{p-k,-m2n,-2p + k,+l)\k,\{nf+p-n,y.{2n,-2p-l)\ 



(23) 



which may be written as 



f \ 



■,2«,-l 



Vi=o 



(24) 



pi 
/=1 



V 



(-1) 



V 



2'-2^+2(2n^+/ + l) 



J J 



Carrying out the summation over i and simplifying the product, (24) becomes 



^0 



which is so far an exact expression, 



(25) 



Using Stirling's formula one excellent working approximation to (25) even 
when p is relatively small is 

b, pP"n/'' [2nf + 2pf [2nf+p^ ^jf^pn.V4 

(eWf(2nf+pf'"''' 
(or often more usefiiUy written in its log form 
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( ( 



exp 



log 



"/■log 



Z?o p^'^n/'^ (lUf + 2pf [in J +p + l) 
(e^p'f^lnf+pf^" 

[luf +2p^{2nf +/7 + 1) 



p/2+1/4 



(iHf+p) 



(26) 



to avoid terms of order Hj"^ ). 

A simpler but more approximate expression can be used for calculations that 
involve large values of and p , provided n, » p»\. Then 

ni ^nf[=ni - p) , and using lim ((2«-/7 + l)/(2«-2/7))" =exp((/7 + l)/2) , 

enables (26) to be simplified to 



{ ^ 1 


4 




2 


f \ 
P. 













£-3 
2 



(27) 



The usefulness of the approximate equations (26) and (27) in determining 
^R~\o ^nfifi^)'' ^niiii^)'^^ tested by comparison with the direct 

integral using lower values for n,- and p and this is shown in table 1 in units of 

Table 1. Value of Ijf^ = Rjjj-ij-(r)r Rjjiif(r) dr obtained by (a) direct 

integration, (b) using equation (26), and (c) using equation (27), for an initial 
state «, , Pi and final state ending on a Pf =^ diagonal of figure 1 . 







Radial 


Value using 


Value using 




p 


overlap 


equation 


equation 




integral 


(26) 


(27) 


10^ 


5 


5.76 


5.80 


6.43 


10^ 


10 


5.58x10"^ 


5.60x10"^ 


7.37x10"^ 
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10^ 


50 


6.88x10"^ 


6.89x10"^ 


6.08x10"^ 


3x10^ 


5 


6.38 


6.43 


6.43 


3x10' 


10 


4.70x10"^ 


4.72x10"^ 


4.73x10"^ 


3xl0' 


50 


1.10x10"^^ 


1.10x10"^*^ 


1.12x10"^^ 



The approximations (26) and (27) highlights the rapidity with which 
Ir-\q ^nf,if^^)^ ^nifii'^) ^'^ approaches zero as p increases, because of the 



factor 



p-3 



2 

. Table 2 provides values of for a range of values of p and «, 



illustrating this effect. The overwhelming smallness of /^j and the resulting 
overwhelming smallness of the dipole matrix element means that state -to- 
state decay rates which start on any large-/) diagonals and end on the p-l 
diagonal (or as it turns out on any low-p diagonal) are negligibly small, even 
despite sizeable values of n^-nj and the cubic dependence of Ajj on . 

This is an extremely interesting result since it means that all the high-n, low-p 
eigenstates in gravitational wells discussed here all exhibit the unusual property 
of having extremely long lifetimes. As a consequence they will, despite their 
high-valued quantum parameters, also exhibit inherent long-term stability and an 
inability to coalesce by radiative decay or interact strongly with other 
electromagnetic radiation. 

Table 2. Value of Ij^ =J^ R^y i^(r)r R^m{r) dr obtained using equation (26), 

for an initial 'deep' state D (quantum parameters , »1) that ends on the 
p- \ diagonal of figure 1. 



\^Kf,ifir)r'R„i^uir)dr 
(in units of Zjq) 

1000 1 -10* 

1000 5 ~6 
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1000 
1000 



10 



30 



10 



30 



10 



30 



10 



30 



10 



30 



X 10^^ 5 X 10 



20 
100 

1 

5 

10 



10 



,20 



10 



26 



.31 



-3 xlO 



-11 



40 

40' 
~6 
40"' 



■40 



60 



40 



-5x10" 



10 



10 



-2x10"' 



-5x10^ 



4. Conclusion 

Exact general expressions have been determined for the transition probability 
per unit time of the electromagnetic decay of charged particles in gravitational 
eigenstates for any dipole-allowed transitions in weak gravity. Additionally 
tractable approximations have been obtained for certain specific transitions 
involving the interesting set of high-n, Xow-p eigenstates, including a general one 
for transitions ending on the p ^ I diagonal. The generality of these types of 
formulae also give them application in the study of other large-« quantum 
systems such as Rydberg states in atomic systems. The most significant 
conclusion to come out of this present study is that the very high-«, low-p 
gravitational eigenstates have extremely long lifetimes (potentially many times 
the age of the universe in the situations considered here). 

The speculation was raised in the introduction that gravitational eigenstates 
might exist naturally elsewhere in the universe. Small structures in low-central- 
mass wells are generally weakly bound (as demonstrated in [1,2]), so that if well 
bound eigenstructures exist they are likely to be of macroscopic size. This is not 
a problem in itself, but it is hard to see how these structures could have formed 
in all but the earliest times in cosmic history since it can be shown [4] that the 
probability for a localised object of any reasonable mass decaying into an 
eigenstate in any time is negligible. The formation of such structures would 
require the existence of strong gravitational potential wells at a time in universal 
history when the global particle density over the well region was sufficiently 
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high and the temperature sufficiently low to enable decay into eigenstates, 
analogous to plasma recombination. Indeed if such conditions were met it would 
be expected that eigenstates would form and ultimately populate the longest 
lived low-/> levels. It turns out that such conditions are possible at certain critical 
times in universal history and the formation mechanisms associated with these 
structures will be the subject of a later paper. It is clear however that the 
resulting stability of particles occupying eigenstates would mean that if such a 
structure formed, it would be stable and have an inability to gravitationally 
collapse. It is also expected that it would not intrinsically emit or scatter 
radiation to any significant degree and therefore be essentially invisible with 
respect to external electromagnetic radiation [4]. This latter point will be 
discussed further in a companion paper. 

The approximation for calculating the dipole transition rate for those transitions 
involving p = \ is a useful one, but it is clearly necessary to develop further 

approximate techniques for the calculation of other dipole transitions with more 
generahzed values of p. In this regard, the derivations in the present paper 
continue to have limitations. In the case of the gravitationally bound systems 

with quantum parameters n, I > 10^" for example, the direct use of (20) is 
limited because of the very large factorial functions (despite the use of Stirling's 
approximation) and also because it involves summations with unreahstically 
large numbers of terms. Different approaches will be needed to deal with these 
types of transitions, where both the initial and final p values of the states are 
large. 

Other areas requiring more detailed examination are those of multi-pole decay 
and transitions induced by other types of interactions such as particle collisions. 
These investigations will be the subject of later papers. It would appear however 
from this initial work that the high-«, low-/) gravitational eigenstates are 
particularly interesting to study theoretically because of their extremely long 
lifetimes and the implications this has for their expected behaviour and 
appearance. 
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